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ABSTRACT 

Aims. We study the characteristics of field stars in the two bright ellipticals NGC4621 and NGC4374 in the Virgo 
cluster to derive distances and stellar population properties. Moreover, since the target galaxies have hosted three type 
la Supernova events, we investigate the correlations between the SNela properties and their host stellar systems. 
Methods. Using deep imaging BVR data, obtained with the FORS2 camera mounted at the VLT, we analyse the 
Surface Brightness Fluctuations (SBF) properties of the targets. We adopt our measurements and existing empirical 
or theoretical calibrations to estimate the distance of the NGC 4621 and NGC 4374. For stellar population analysis, 
we measured SBF amplitudes in different galaxy regions, allowing to study the changes in field star properties. The 
three band coverage of present data, coupled with existing SBF measures available from literature, provides us with the 
largest wavelength coverage of SBF magnitudes for single objects. We present a detailed comparison between SBF data 
and models to constrain the physical characteristics of the dominant stellar components at i) various galactic radii, and 
ii) in the regions where SNela events were recorded. 

Results. Our V and R SBF measures provide distances in agreement with literature estimates. The median of our and 
literature SBF-based distances agrees with the one from non-SBF methods, indicating the absence of any systematic 
effect in the SBF technique. Comparing either the SBF versus integrated colour diagrams, or the SBF versus SBF 
colour diagrams, with models, we find that stellar populations properties do not change significantly along galactic 
radius, with a dominant population having old age and solar chemical composition. The galaxies appear similar in all 
properties analysed, except for B-band SBF. Since the SBF magnitudes in this band are sensitive to the properties of 
the hot stellar component, we speculate that such behaviour is a consequence of different diffuse hot stellar components 
in the galaxies. Using specific models we find that the presence of a percentage of hot-HB stars in old and metal rich 
stellar populations could be at the origin of the observed differences. We find a good uniformity in the V and R SBF 
and integrated colours in the regions where the three SNela, presenting different absolute luminosities, exploded. On 
the other hand, the B-band SBF signal shows intriguing differences. 

Key words. Galaxies: distances and redshifts - Galaxies: elliptical and lenticular, cD - Galaxies: stellar content - 
Galaxies: photometry - Galaxies: individual: NGC 4621, NGC 4374 



. !_h ■ 1. Introduction galaxies, has proved being not only an accurate and precise 

distance indicator, but also a powerful tracer of stellar pop- 

; Our understanding of galaxies m the distant Universe re- u i atio n properties fe.g. I Jensen et al. 2003; R aimondo et all 

Ctf . lies on how well we understand the properties of their lo- gj7j| ICantiello et al.ll2007bi . Bv definition, the SBF mag- 

cal templates. Thus, the study of nearby unresolved stellar nitude corresponds to the ratio of the second to the first 

populations p ays a key role to obtain a refined character- moment of the stellar i uminosi ty function in the galaxy: 

isation of stellar populations at larger redshifts. Although _ - s^/ 2 

different astronomical methods exist to carry out such de- m = -2.5 log/, where / = f. , and m is the n umber 

tailed analysis, none of them can provide by itself robust of stars per bin of flux /j (jTonrv k Schneiderl [ l988h . Such 

constraints. The presence of internal uncertainties in each definition implies that: (i) SBF magnitudes are linked to 

method, or calibration uncertainties, as well as the effect the properties of the stars in the galaxy; (ii) the SBF sig- 

of the age-metallicity degenera cy that affects most of the nal is dominated by the brightest stellar component in the 

spectro-photometric indicators (jWorthevll 19941 ). prevents us galaxy, because of the dependence on the second moment 

from relying on one single stellar population tracer for these of the luminosity function; (Hi) since the brightest phase in 

studies, and pushes the community of astronomers towards a stellar population is wavelength-dependent, SBF magni- 

the definition of new analysis methods. tudes in different pass-bands are sensitive to the properties 

In the last two decades , the SBF method, introduced by of stars in different, and well defined, evolutionary stage s 

iTonrv fc Schneider! (|1988T ) to obtain distances of early-type (|Worthevlll993al iBrocato et aTlll998t ICantiello et aD2003[) . 
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Cantiello et al.: BVR SBF analysis for NGC4621 and NGC4374 
Table 1. Main parameters of the target galaxies and observations 





NGC4621 


NGC4374 


Galaxy parameters 


RA(J2000) i 
Dec(J2000) 1 


12h42m02.3s 


12h25m03.7s 


+ Ild38m49s 


+12d53ml3s 


Galaxy Type 2 E5 
Morphological Type 2 -4.8 
Absolute B-band magnitude 2 -20.5 


E 
-4.3 
-21.0 


SNe la events 


SN1939B 


SN1957B, SN1991bg 


cz 1 (km/s, Heliocentric) 

E(B-Vf 

{V-If 


410±6 
0.033 mag 
1.172 ± 0.018 


1060 ±6 
0.040 mag 
1.191±0.008 




Observations 




Filter 
B 


Exposure time (s) 
3375 


Exposure time (s) 
2250 


V 


3375 


2250 


R 


2376 


1584 



1 Data retrieved from NED (http://nedwww.ipac.caltech.edu); 2 Hyperleda (http://leda.univ- lyonl.fr [ ); 3 ISchleeel et~aD (| 199S 
4 iTonrv et~aTl (|200ll ) 



Taking advantage of archival B, V , and R observations 
taken with the FORS2 camera of the Very Large Telescope 
for two bright early- type galaxies in the Virgo cluster, 
NGC 4621 and NGC 4374, we obtain multi-band SBF mea- 
surements of these galaxies. These measurements, coupled 
with existing ground-based I— and K— band measures, and 
space-based F850LP from ACS data (z hereafter; for inte- 
grated colours we will also use ACS F475W data, indicated 
as g in the paper), provide the largest wavelength coverage 
of SBF measures for single galaxies. With such sample of 
measures we investigate the properties of the two galax- 
ies using the SBF technique under its twofold aspects: as 
a distance indicator and as a stellar population properties 
tracer. Furthermore, since the target galaxi es have hoste d 
three SNe la eve nts - SN 1939B in NGC 4621 (|Zwicky|1942fl , 
and SN 195 7B dBertola 11964 . and references therein) plus 
SN1991bg (iKosai et al.lll991h in NGC 4374 - we can ex- 
plore the capabilities of the SBF method to improve our 
knowledge on the relation between the SNe la progenitors 
and the stellar population they belong to. 

The paper is organised as follows: a description of the 
observational data set, the data reduction and calibration 
procedures, and the procedure for SBF measurements is 
given in fj3J The analysis of distances is presented in 521 
while the study of stellar population properties based on 
SBF is presented in 2) A summary and the conclusions 
finish the paper in [j5l 



2. The Data 

For this work we used B, V and R data of two Virgo clus- 
ter galaxies, NGC 4621 and NGC 4374, observed with the 
FORS2 camera of the VLT. Some relevant properties of the 
targets are reported in Table [1] 

The science images of the galaxies and the calibration 
files (bias, fiat, field of standard stars) were retrieved from 
the ESO archivefl 

In one case, NGC 4374, we could not use the total inte- 
gration time available because of a slight rotation ( ^ 1 deg) 



between different frames which badly affects the SBF signal 
(see below). In Table Q] we list the total exposure times for 
each filter and for both galaxies. 

2.1. Data reduction and calibration 

All data were retrieved from the archive and processed with 
standard data reduction procedures using IRAF0 tasks. 
The procedure is briefly described here. We obtained mas- 
ter bias and flat frames (one per filter) for each observ- 
ing night. Individual frames are bias-subtracted, flattened, 
shifted into registration, and combined, rejecting pixels 
affected by cosmic ray hits. To combine the frames we 
imposed no sub-pixel registration, in order to avoid the 
contamination to S BF due to the sub-pixel interpolation 
(jJensen et al.lll998h . Finally, the seeing between targets and 
filters ranged between 0".8 and 1".0. The combined frames 
are shown in Figure [1] 

The standard calibration plan of FORS2 pr ovides 
nightly multi-band observations of iLandold (|1992[ ) stan- 
dards, which are used to calibrate the photometry of the 
two galaxies. 

2.2. Data analysis and SBF measurements 

To derive the photometry of sources in the frames and mea- 
sure the fluctuation amplitudes we used the same proce- 
dure d esc ribed i n our previous works (jCantiello et al.l l2005. 
l2007allbl 120091: iBiscardi et ail I2008D . The main steps of 
SBF measurement involve: sky background determination, 
model and large scale residual subtraction; photometry and 
masking of point-like and extended sources, including dust; 
power spectrum analysis of the residual frame. 

Here we briefly describe some relevant parts of the anal- 
ysis. We determined the sky background in the galaxy im- 
ages by fitting the surface brightness profile of the galaxy 



http:/ /archive. eso.org/ 



2 IRAF is distributed by the National Optical Astronomy 
Observatories, which are operated by the Association of 
Universities for Research in Astronomy, Inc., under cooperative 
agreement with the National Science Foundation. 
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r r 

NGC4621 NGC4374 

Fig. 1. Left Panel: NGC4621 ^-band image with the regions used for SBF and colour analysis. Thick-lined circles mark 
the borders of the main annulus, thin solid lines show the other annuli (as numbered), white boxes give the location of 
the SNe (see £|2.2j) . Right panel: as upper, but for NGC4374 (for sake of clarity the inner border of the main annulus is 
shown with white colour). 



with a Sersic's law (jSersid [l968h plus a constant term. A 
first model of the galaxy was obtained and subtracted from 
the sky-subtracted frame. The wealth of bright sources ap- 
pearing after model subtraction were masked out; the pro- 
cedure of model fitting and masking was then iteratively re- 
peated until the residual frame (original frame minus galaxy 
model) was considered satisfactory. After the subtraction of 
the best galaxy model, the large scale residuals still present 
in the frame were remo ved using the backgroun d map ob- 
tained with SExtractor (jBertin fc Arnouta 1996) adoptin g 



a mesh size ~ 10 times the FWHM (|Cantiello et al. 2005). 
In the following we refer to the sky, galaxy-model and large 
scale residuals subtracted image as residual frame. 

The photometry of fore/background sources and of 
Globular Clusters (GCs hereafter) was derived with 
SExtractor on the final residual frame. As described in our 
previous works, we modified the input weighting image of 
SExtractor by adding the galaxy model (times a factor be- 
tween 0.5 an d 10, depending on the galaxy; for details see 
Uordan|[200l ICantiello et al.ll2005b so that the SBFs were 
not detected as real objects. To correc t for Galactic redden- 
ing we used the E(B—V) values from lSchlegel et all (|1998f ) 
reported in Table [TJ The aperture correction (a.c.) was ob- 
tained from a number of isolated point-source candidates 
in the frames and by making a curve growth analysis out 
to 6.0". 

Once the corrected catalogue of sources was derived, 
the next step was to fit the luminosity functions of external 
sources, that will be used to estima te the extra-fluctu ation 
term due to unmasked faint sources (|Tonrv et al. 199C|). We 
derived the fit to GCs and background galaxies luminosity 
functions from the photometric catalogue of sources, after 
removing the brightest / saturated point-like sources and the 
brightest and most extended objects. The best fit of the 
sum of the two luminosity functions, shown in Figure [21 




20 22 24 26 20 22 24 26 20 22 24 26 
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Fig. 2. B, V, and R luminosity functions of external sources 
in NGC4621 (upper panels), and NGC4374 (lower panels). 
Open squares mark observational data; the solid lines repre- 
sent the best fit to data. The luminosity functions of back- 
ground galaxies and GCs are plotted as dotted and dashed 
lines, respectively. 



and t he extra-fluctuation c orrection term, P r , were derived 
as in ICantiello et all (l2005l) . 

To measure SBF magnitudes we proceeded by estimat- 
ing the azimuthal average of the residual frame power 
spectrum, P{k), then matching it with the power spec- 
trum of a template PSF convolved with the mask image, 
E(k). The total fluctuation amplit ude Pg was obtai ned 
via a robust minimisation method (jPress et al.l Il992l ) as 
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the multiplicative factor required to obtain the matching 
P(k) = Pq x E(k) + Pi, where Pi is the white noise con- 
stant term. As template PSF we used 6 different isolated 
bright point-like sources in each residual frame, which, af- 
ter normalisation, were singularly adopted to estimate the 
SBF signal of the galaxy. 

For both objects the fluctuation amplitude, Pf = Pq — 
P r , was estimated within a circular annulus, adopting the 
same average inner and outer radii used by iTonrv et al.l 
(200] ]) in order t o allow a homogeneous combination of our 
and ITonrv et al.l data. Figure [3] shows the power spectrum 
analysis of both galaxies. 

Furthermore, in order to study the radial variation of 
SBFs, and the SBF properties in the sites where Type 
la Supernova events are recorded, we ran the same SBF 
measurement procedure described above in five annuli per 
galaxy, and in three more box-shaped regions centred on 
each SN. All regions used are shown in Figure [1] We must 
emphasise that the effect of the extra-fluctuation correction 
term, P r , and its relative amplitude with respect to the to- 
tal fluctuation signal, Pq, changes from region to region, and 
the P r /Po ratio increases for outer regions and bluer bands. 
Thus, a systematic under/over-estimate of P r may lead to 
over /under- estimated SBF magnitudes. Table [5] shows the 
(Pr/Po) x 100 ratios for all regions and bands considered 
here; the numbers in the table demonstrate that V and R- 
band SBF magnitudes have P r corrections much lower than 
the total fluctuation signal, even for the outermost annu- 
lus considered. On the contrary, for B-band measures the 
P r /Po ratio in some of the regions considered reaches val- 
ues as high as ~ 0.45, implying that the extra correction 
term is comparable with Pf, thus the SBF magnitudes will 
be considered with caution in such regions. 

Table [3] summarises our measurements; for each galaxy, 
we tabulate (1) region identification; (2) average equiva- 
lent radius of the region r(arcsec) = y/(a ■ b); (3-5) SBF 
magnitudes measured in BVR; (6-7) integrated colours. 



3. Analysis of distances 

To estimate distances based on the SBF method, as for 
most distance indicators, a calibration of the absolute SBF 
magnitude is required. The most widely adopted bands for 
SBF measures are the optical I and z bands, and the near- 
IR H and K ones dPahre fc Mouldl|l994tlTonrv et al.ll200lt 
iMei et al.l [2007t Uensen et all 12003: IBlakeslee et all 120101) . 
All z and H measures available come from HST observa- 
tions taken with the A CS/WFC and NICMQS/NIC2 detec- 
tors, respectively (e. g. Uensen et all 120031: IMei et all [20071: 
IBlakeslee et~aT1l2009D . The total number of SBF data avail- 
able in these pass-bands is around ~ 450 measures. As a 
consequence, the calibrations in such bands, especially the 
optical ones, are well established as testified, for example, 
by the accurate characterisation of the z versus (g—z) re- 
lation derived by IBlakeslee et al.l (|2009f ) using ACS data of 
Virgo and Fornax cluster galaxies. 

SBF measurements in other bands are not popular, with 
a total of ~7 measures available for B V and R bands (e.. 



less reliable, while for distance studies the most favourable 
condition is that the magnitude of the standard candle 
is relatively constant or has a tight predictive relation as 
a function of some other distance— independent property, 
such as colour. Taking advantage of the well constrained 
distances of the target galaxies, we will derive new distances 
based on our measurements and empirical/theoretical cal- 
ibrations, and compare the various results to test calibra- 
tions adopted. In the upper part of Table 2] we report the 
SBF measures for NGC4621 and NGC4374 in different 
bands as derived in literature, the distance moduli m — M, 
and the calibration used. The second part of the table gives 
the distances derived by us as discussed below. 

3.1. Empirical calibrations 



a total or ~ tu measures available lor tt V ana K bands (e.g . 
Tonrv et al.lll990l: Uerien et alj|2000t IBlakeslee et al.l l200jj: 



Sodemann fc Thomsenl 119961: ICantiello et al.l l2007bh . This 
is due to the fact that SBF magnitudes at these wavelengths 
tend to be more sensitive to stellar population properties 
and, consequently, the calibration has a larger scatter and is 



ICantiello et al.l (|2007aD compared the existing / and V 
band calibration equations derived from different observa- 
tional data, in order to identify the best empirical calibra- 
tion for the colour interval 0.85 < V — I < 1.30 mag in both 
pass-bands. For F-band SBF magnitudes the best calibra- 
tion equations resulted to be: 

My = 0.81 ± 0.12 + (5.3 ± 0.8) • [(V-I) - 1.15] for 1.05 < 

V-I < 1.30 mag (1); 

and 

My = -0.50 ± 0.27 mag for the 0.85 < V-I < 1.05 mag 

(2) . 

T o use the above cali bration we adopt the V—I colours 
from ITonrv et al.l (|200lD . and the SBF measures listed in 
Table [3] for the main annulus. The distances obtained are 
shown in Table [4] and will be discussed with all other esti- 
mates at the end of the section. 

For i?-band SBF, the most recent empiric al calibration 
for no rmal elliptical galaxies was provided bv ITonrv et al.l 
(1990) who, however, considered it unreliable. Distance esti- 
mates for a large number of dwarf ellipti cals have been pro- 
vided by Je rjen and collaborators (e.g. iJerien et al.l il998. 
120001 l200l based on i?-band SBF measures and semi- 
empirical relations. Using the B—R colour, the authors 
recognised two different branches for the SBF versus colour 
relation, a linear branch: 

M R = -8.94 + 6.09 • (B-R) for 1.10 < B-R < 1.50 mag 

(3) ; 

and an parabolic branch, partly overlapping to the previous 
one: 

M R = - 1.39+1.89- [(B-R) -0.77] 2 for 0.80 < B-R < 1.35 
mag (4) (perien et al.ll2004[) . To derive distances, reported 
in Table HI we used the red/linear branch calibration. 

The small amplitude of the fluctuations in the B-band 
and the higher sensitivity to stellar population properties 
with respect to other opti c al ban ds (see, e.g., IWorthevI 
Il993at ICantiello et all 120031. l2007bh make this band unre- 
liable for distances studies. Besides, only a handful of B 
measures exist, including the present ones, and no distance 
has been estimated using B-band SBF measurements. 

3.2. Theoretical calibrations 

Empirical calibrations are generally preferred for SBF stud- 
ies, also thanks to the aforementioned small uncertain- 
ties for certain filters. However, the derivation of such re- 
lations requires substantial observational effort for each 
pass-band. On the other hand, models have the advan- 
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Fig. 3. Power spectrum analysis of the main annuli for both targets. Left panels: NGC4621 analysis. The upper panel 
shows the logarithm of the power spectrum of the residual frame (grey dots) and the best fit P(k) curve. In middle panel 
the difference between observed and fitted data is shown. The flat region of log Ppjkfjt > k) between vertical dashe d 
lines in the lower panel is used to evaluate the best fit parameters Po, and Pi (jCantiello et al.ll2005t iBiscardi et al.l l2008). 
Right panels: as left, but for NGC 4374. 

Table 2. (P r /P ) ratios (in %) 





B 


V 


R 


B 


V 


R 


Anmilus 


NGC 4621 


NGC 4374 


Main 


17 


7 


5 


29 


11 


9 


1 


6 


4 


4 


11 


7 


8 


2 


13 


6 


4 


17 


9 


8 


3 


24 


9 


6 


26 


11 


9 


4 


34 


13 


8 


38 


13 


9 


5 


37 


15 


10 


44 


16 


11 


SN 1939B 


12 


5 


4 




SN 1991bg 








37 


10 


9 


SN 1957B 








25 


9 


7 



Table 3. Surface Brightness Fluctuation and colour measurements corrected for galactic extinction. 



Annulus 


(r> 


B 


V 


R 


B-V 


B-R 




(arcsec) 


(mag) 


(mag) 


(mag) 


(mag) 


(mag) 








NGC 4621 






Main 


81 


32.83 ± 0.10 


31.84 ± 0.06 


30.94 ± 0.08 


0.893 ± 0.002 


1.449 ± 0.002 


1 


29 


32.43 ± 0.09 


31.75 ± 0.07 


30.96 ± 0.06 


0.912 ± 0.001 


1.474 ± 0.001 


2 


58 


32.83 ± 0.10 


31.83 ± 0.06 


30.94 ± 0.08 


0.900 ± 0.001 


1.455 ± 0.001 


3 


95 


32.83 ± 0.10 


31.81 ± 0.06 


30.93 ± 0.07 


0.884 ± 0.002 


1.439 ± 0.003 


4 


132 


32.87 ± 0.11 


31.74 ± 0.07 


30.90 ± 0.08 


0.871 ± 0.005 


1.431 ± 0.006 


5 


193 


32.63 ± 0.12 


31.59 ± 0.06 


30.81 ± 0.06 


0.839 ± 0.008 


1.391 ± 0.011 


SN 1939B 


55 


33.01 ± 0.10 


31.84 ± 0.07 


30.97 ± 0.08 


0.915 ± 0.001 


1.469 ± 0.001 








NGC 4374 






Main 


93 


33.31 ±0.08 


32.11 ± 0.04 


31.40 ±0.08 


0.925 ± 0.002 


1.494 ± 0.003 


1 


20 


32.43 ±0.06 


31.51 ± 0.04 


31.01 ±0.08 


0.945 ± 0.001 


1.521 ± 0.001 


2 


39 


33.08 ±0.07 


32.03 ± 0.04 


31.33 ±0.08 


0.931 ± 0.001 


1.501 ± 0.002 


3 


63 


33.36 ±0.09 


32.23 ± 0.03 


31.45 ±0.08 


0.922 ± 0.002 


1.489 ± 0.004 


4 


88 


33.58 ±0.08 


32.23 ± 0.02 


31.41 ±0.08 


0.919 ± 0.003 


1.486 ± 0.006 


5 


122 


33.47 ±0.10 


32.22 ± 0.03 


31.44 ±0.08 


0.919 ± 0.006 


1.486 ± 0.011 


SN 1991bg 


58 


33.87 ±0.07 


32.13 ± 0.04 


31.43 ±0.08 


0.927 ± 0.002 


1.495 ± 0.00 


SN 1957B 


59 


33.50 ±0.09 


32.13 ± 0.03 


31.35 ±0.08 


0.917 ± 0.002 


1.484 ± 0.00 
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tage of being homogeneous through the different bands, 
but need many different counter-checks to be consid- 
ered reliable. In this work we take as reference the SBF 
versus colour equations derived using the simple stel- 
lar population (SSP) models from the Teramo-Stellar 
Population Tools (SPoT group0. For a d etailed review of 
the SP oT models we remind the reader to iRaimondo et al.l 
(2005) and references therein. These models have already 
been proved being very effective not only to match the 
empirical SBF calibrations in different bands, but also 
to reproduce the resolved (colour magnitude diagrams) 
and unresolved (colours, magnitudes) properties of stel- 
lar populations jlBrocato et alJl200(it [Cantiello et aL 2003: 
Raimondo ct al. 2005; Cantiello et al Jl2007at iBiscardi et aJj 



Keeping in mind such uncertainty, and the error in the 
SBF measures, the results in Table [4] can be summarised as 
follows: 



2008t ICantiello et al.ll2009t lRaimondoll2009D 



For example, in a previous paper. [Biscardi et al.l (|2008| ) 
showed that the SPoT models for an age interval between 
1.5 and 14 Gyr, and metallicity [Fe/H] between -0.4 and 
0.3 dex are able to reproduce both t he J -band and z-ban d 
calibrations from lTonrv et al.l (|200ll) and lMei et all (j2007lh 
respectively. 

In the present work we adopt the same models used 
by IBiscardi et al.l and derive absolute V and i?-band SBF 
magnitudes versus the integrated B—R colour. Using the 
on-line SPoT models we obtain the following relations: 
M R = -4.11 + 2.84 • {B-R) (5), 
My = -2.69 + 2.60 • (B-R) (6), 

we also derive, for sake of completeness, the -B-band cali- 
bration: 

M B = -0.68 + 2.32 • (B-R) (7). However, we recall that_it 
is hazardous to derive distances from eq. (7) because Mb 
strongly depends on the properties of the stellar population 
originating the SBF signal (see section § 14. 1|) . 

As a further check based on an independent set of stel- 
lar populati on models, we also use d the V and R equations 
provided bv lBlakeslee et ail (|200lD . which, differently from 
the SPoT models, are derived using composite stellar pop- 
ulations (see the quoted paper for more details). 

The distance moduli of the two targets, ob tained using 
both the SPoT and the iBlakeslee etall (|200lD theoretical 
calibrations, are reported in Table 0J 

3.3. Results 

By inspecting the new and old distance estimates in Table 
0] we find a satisfactory agreement within the quoted un- 
certainties, no matter what calibration (empirical or theo- 
retical) or pass-band/colour relation is used, with the sole 
exception of _B-band data which we report here only for 
sake of completeness. 

The distance moduli derived from SBF are affected by 
the uncertainties present in the empirical/theoretical cali- 
brations, besides the uncertainty of the SBF measure itself. 
In this work, for the distance moduli taken from literature 
where no calibration error is given we consider safe to as- 
sume an error of the order of ~0.20 mag, which includes 
zero-p oint uncertainty and the scatter of em pirical calibra- 
tions (|Tonrv et al.l 120011 1 Jensen et al.ll2003l ). Similarly we 
assume ^0.20 mag error for the theoretical calibrations, 
which is originated by the spread between models with dif- 
ferent age and metallicity. 



— although the colour provided by iTonrv et al.l was mea- 
sured in galaxy regions slightly different from our fgthe 
distances derived from V using eq. (1) agree with other 
data from literature. A simila r result is tr u e for R if 
the linear branch relations by iJerien et al.l (|2004f ) are 
adopted; 

— the distance estimates obtained using empirical calibra- 
tions show a slight tendency to have a smaller scatter 
with respect to those from theoretical calibrations; 

— whether SSP models from the SPoT group or the 
IBlakeslee et al.l (|200l|) composite models are used the 
distances derived are substantially similar to each other, 
and agree with expected values. Such result implies that 
for the bandpasses and the colour interval considered 
here composite stellar population models, which try to 
better reflect the real population mixing of galaxies, are 
not strictly necessary for the purposes of deriving ap- 
propriate SBF versus colour relations; 

— as discussed above, coupling B measures with eq. (7) 
provides unreliable distance moduli, this confirms that 
this band must discarded for distances. Interestingly, 
the difference between the i3-band distance moduli and 
the literature average moduli is significantly larger for 
NGC4621 (A(m-M) = (m-M) B ,SBF-(m-M) Ht . ~ 

-0.8 mag) than for NGC4374 (A(m-M) 0.6 mag). 

We will analyse further this evidence in the next section 
on stellar population analysis. 

Finally, the last lines of Table 2] provide the median dis- 
tance moduli of the two targets obtained by averaging i) 
our estimates (empirical and theoretical calibrations, ex- 
cept i?-band), it) all SBF m — M, including ours, Hi) the 
distance moduli without SBF^I, and iv) the results from 
ours plus all literature data. By inspecting such data we 
find an excellent agreement between our and literature dis- 
tance moduli and, more in general, between non-SBF and 
SBF-based distances, a results that should be regarded as 
a direct evidence of the absence of any significant bias or 
systematics between the quoted distances. In addition, such 
agreement also proves the reliability of V-and i?-band SBF 
calibrations presented in this section. 

3.4. Comparison with SNela distances 

The comparison of SBF and SNela distances plays a fun- 
damental role in the problem of the cosmological distance 
scale. SBF magnitudes are capable of providing distances 
with < 10% accuracy within 100 Mpc dJensen et al.ll200lt 
IBiscardi et al . 2008; Blakc slee et al.ll2010l ). a result that will 



Visit the web-site: www.oa-teramo.inaf.it/spot 



4 As explained b efore, we ado pted the same "average" inner 
and outer radii of ITonrv et al.1 . but the detailed shape of the 
annulus, plus the masking of sources is clearly different between 
the two datasets. 

5 All non-SBF distances are taken from the NED Redshift 
Independent Distance database. The distance indicators used 
include Globular Cluster and Planetary Nebulae Luminosity 
Functions, Type la Supernovae and Globular Cluster half light 
radii. For a complete list of methods, and associated references, 
visit the URL http://ned.ipac.caltech.edu/forms/d.html. In case 
of multiple estimates obtained with the same indicator only the 
most recent one is considered. 
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Table 4. Galaxy distances from this work and from literature. 











NGC4621 


NGC4374 


Filter 


Colour Ref. Color/Data 


m m — M 


m 


m-M 


SBF from literature 


I 


V 


-I 


W 


29.67± 0.18 31.16±0.20 


29.77± 0.09 


31.16±0.11 


z 


g- 


-z 


[2] 


29.12 a ± 0.01 30.86±0.07 


29.53 a ± 0.01 


31.34±0.07 


K 






[3] 


25.46± 0.16 


25.43± 0.22 




Our measurements 


Empirical Calibrations' 3 


V 


V 


-I 


[4j 


31.84± 0.06 30.91±0.21 


32.11± 0.04 


31.08±0.20 


R 


B- 


-R 


[5] 


30.94± 0.08 31.05±0.21 


31.40± 0.08 


31.25±0.22 


Theoretical Calibrations-SSP b 


B c 


B- 


-R 


L6J 


32.83 ±0.10 30.2 


33.31 ±0.08 


30.5 


V 


B- 


-R 


[6] 


31.84 ±0.06 30.76±0.21 


32.11 ±0.04 


30.91±0.20 


R 


B- 


-R 


[6] 


30.94 ±0.08 30.93±0.21 


31.40 ±0.08 


31.27±0.22 


Theoretical Calibrations-CSP b 


V 


V 


-I 


m 


31.84± 0.06 30.99±0.21 


32.11± 0.04 


31.17±0.20 


R 


V 


-I 


[7] 


30.94± 0.08 31.04±0.21 


31.40± 0.08 


31.41±0.22 


Median m — M 










NGC4621 


NGC4374 






This work 




30.96 ±0.11 


31.21± 0.17 






all SBFs 




31.08 ±0.19 


31.21± 0.14 






No SBFs 




30.97 ±0.21 


31.17± 0.18 




All available m — 


M 


30.98 ±0.17 


31.22± 0.16 
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[1 ] iTonrv et alJ d2001f>: [21 IBlakeslee et all (|2009h : [3] iPahre fc Mouldl l| 19941 ): [4] ICantiello etall (|2007al 1: [5] Uerien et alJ lj2004Tl : 
[6] this work; [7] IBlakeslee et all (|200ir K 

Notes 

a IBlakeslee et alJ (|2009T ) data are given in the AB-mag system. 

b We adopted ~0.2 mag as safe calibration error for i) distance moduli given in literature with no calibration uncertainty, and 

ii) distances based on theoretical calibrations (see text), 
c Distances derived from B-band SBF are reported to show that the sensitivity of B to stellar population poperties makes 

unreliable any distance estimate. 



be likely improved with new and next generation observing 
facilities, allowing to estimate accurate distances of bright 
ellipticals out to ^200 Mpc. On the other hand, SNela can 
provide distances to objects at much larger distances. Thus, 
deriving self-consistent distances using these two indicators 
is a crucial step to bridge local to cosmological distances, 
in order to reduce the number of rungs in the cosmological 
distance scale, i.e. to reduce systematic uncertainties. 

The sample of objects with known SBF measures and 
well studied SNela light curves is relatively small, due to 
the fact that SNela occur in all kind of galaxies, but are 
mostly observed in late type galaxies because of an observa- 
tional bias, on the contrary SBF are measured almost only 
in early type galaxies. 

Nevertheless, a c omparison of S BF an d SNela distances 
was carried out by lAihar et al.l (|200lt ). who found that 
there is a good agreement between the two distance in- 
dicators, provided that a consistent set of Cepheid cali- 
brating galaxies is us ed. However, their statistics grounded 
on ~ 10 objects and lAihar et al.l made it clear that "[...] 
Unquestionably, the SNIa and SBF absolute calibrations 
are in need of further refinement. " 

As already mentioned, a total of three type la Supernova 
events has been recorded in the two galaxies. Two SNIa 
have been observed in NGC 4374, SN 1991bg and SN 1957B, 
while SN1939B was discovered in NGC 4621. The SNe are 



located all at ~1 arcmin from the photometric centre of the 
host galaxy. Even though all three SNe are classified as type 
la, they show a w ide range of luminosit y - not unexpected 
in E/S0 galaxies (Gallagher et al. 2005) -, in particular all 
events are fainter at maximum light than a normal SNIa 
after correcting for t he correlation betw een peak luminos- 
ity and decline rate (jHamuv et al.|[T996l ). SN 1991bg is one 
of the most sub- lu minous SNIa yet obs erved (~ 2.5 mag 
fainter than normal. [Turatto et al.|[l996Tl . SN1957B has an 
absolute magnitude at maximum ~ 0.2 mag faint er than 
a no rmal SN, but still brighter than SN 1991bg (jHowelll 
120011 ). Finally SN 1939B at ma ximum light is ~ 0.6 mag 
sub- luminous (Minkowski EM^). 

Unfortunately, being sub-luminous events, the standard 
peak luminosity versus decline rate relation d oes not pro- 
vide g ood distance estimates with these SNe. lAihar et al.l 
(|2001l ). in fact, did not take into account SN 1939B 
and SN1957B, while SN1991bg appeared in their list of 
Supernovae but it was not used for SBF-SNela compari- 
son, it is was instead included "for completeness and for 
future studies of SNela luminosities." . 

More recent studies, however, have provided new meth- 
ods that can be adopted to derive distances with these 
sub-luminous S Ne: the calibra ti on ba sed on decline rate 
and colours by iFolatelli et all (|2010l ). the MLCS2k2 by 
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Uha et all (120071) . and the AC12 method by IWang et al] 
(1200611 . 

Using eq. (6) of iFolatelli et all (|2010h with the best- 
observed fit parameters, we obtain the following distance 
moduli: 



SN 1939B in NGC 4621: taking the peak luminosity from 
the Asiago Supernova Catalogue (|Barbon et al.l [l989'Fl 
m-B^ax = 12.3, a decline rate of Ami n ~ 1.75 estimated 
from the light curve report ed in iLeibundgut et al.l 
(|1991| ), and using eq. (3) from IFolatelli et all (|2010D to 
get the (B max - V max ) colour, we obtain m-M = 31.0, 
with an uncertainty of the order of ~0.5 mag if errors 
of ^0.2 mag, 0.25 days and 0.15 mag are adopted for 
m B ,max, Atoi 5 and (B max - V max ) , respectively; 
SN1957BinN GC4374: ado pting m R . max = 12.20±0.14 
(|Lanor 



b3 ll998l) . Amu 



= 31.1, with uncertainty 
0.1 on both the decline 



1.3 (Howell 2 00l, and (BJ 1 
V max ) Q = 0.01 (using eq. (3) from IFolatelli et al.ll27)10h 
the distance modulus is to — M 
^0.3 mag adopting an error of 
rate and on (B max - V max ) ; 

SN1991bg in NGC 437 4: with ro B .m» f = 14.75 and 
m Vimax = 1 3.95 ± 0.02 dTuratto et al.l ll996l). Atoi 5 = 

19991 ) we derive to — M = 



1.93 ± 0.10 dPhillips et al 
31.1 ±0.3: 



In all cas es the distance modu li derived using the cali- 
bration from IFolatelli et all (|2010fl agree within uncertain- 
ties with the median values reported in Table 2] 

Uha et all (l2007f) develop ed an updated version of the 
MLCS method dRiess et al J 119981) called MLCS2k2, which 
includes new procedures for the if-correction and for in- 
ternal extinction c orrections. Using the light curve param- 
eters of SN 1991bgUlraet_alJ derive m-M = 31.42 ± 0.10 
adopting H — 65 Km s" 1 Mpc -1 . Although such distance 
agrees with the values reported in Table HI the agreement 
becomes certainly better if H = 73 Km s" 1 Mpc^ 1 is 
taken - a value consis tent with the ones t ypically derived 
from SBF distances ()Tonrv et al.l 120001: iFreedman et al.l 
1200 ID - as suggested by the authors; in that case, in fact, 
one has m-M = 31.17 ± 0.10 mag. 

In the case of IWang et al.l (|2006l ) the calibration pa- 
rameter adopted is the AC12, i.e. the B—V colour at 12 
days past optical maximum. The authors find to — M = 
30.35 ± 0.19, based on the properties of SN 1991bg and us- 
ing Hq — 72 Km s _1 Mpc . Such m—M is more than 2-cr 
out the median values in Table|?]obtained using many inde- 
pendent distance indicat ors, and it also disag ree s with the 
values obtained using the lFolatelli et all (120101) or Una et al.l 
(|2007|) calibrations, thus we flag the Wang et al'l estimate 
as unreliable for the case of SN 1991bg. 

In conclusion, we find that the distances of NGC 4621 
and NGC 4374 derived via the light curve properties of their 
SN la agree perfectly with the SBF distances - and with the 
most recent es t imate s from literature - if the cali bration by 
IFolatelli et all (|2010L useful for all three SNe), or Uha et al.l 
(2003, for the case of SN 1991bg) are used. 



4. Analysis of unresolved stellar populations with 
SBF and integrated colours 

As mentioned before, the SBF magnitude of a stellar popu- 
lation corresponds to the ratio of the second to the first 
moment of its luminosity function. Since its first appli- 
cations, such characteristic has suggested the use of the 
SBF method as a powerful diagnostics of stellar population 
properties. The earliest theore tical studies on this subject 
(|Buzzonilll99l: lWorthevllT993allbh were focused on the rela- 
tionship between SBF and age/metallicity, but also showed 
how the fluctuation amplitude in bluer pass-bands, like B, 
could be used as a useful tracer of the hot stellar compo- 
nents in the galaxy. Despite several authors have confirmed 



and further extended, such early findin gs (Blake slee et al, 
120011 : ICantiello et al.ll200l iGonzalez et al.ll2004l : lRaimondo 



' The updated catalogue is available at the URL 
http://graspa.oapd.inaf.it 



2009), a comprehensive and homogeneous analysis of multi- 
band SBF measurements for a large sample of galaxies is 
still missing. The various set of models have clearly shown 
that optical to near-Ii? SBF colours can be very effective 
to lift the age/metallici ty degeneracy w hich affects "classi- 
cal" integrated colours (|Worthevlll994D. Observationa lly, a 
few I — H data exist in literature (| Jensen et al.l l2003). but 
they mostly refer to spatially non homogeneous regions. 
On the other hand, V — I colours have also been studied 
for various galaxies, but they suffer for an age/metallicity 
degeneracy similar to classical colours, although the in- 
formation obtained using such measures is independent 
and complementa ry to the integrat e d colours derived for 
the same targets jTonr v et aljfl990l: iBlakeslee et al.l 120011: 
ICantiello et al.ll2007allbh ~ 

Concerning the use of SBF to study stellar popula- 
tions, using high quality data has allowed to measure SBF 
variations in optical bands within different regions of the 
galaxy. As a consequence, even if V and I SBF magni- 
tudes and colours are not as much effective in removing 
the age/metallicity degeneracy (see below), the measure of 
SBF radial variations has provided useful results on how the 
mean properties of t he dominant stellar p opulation change 
with galactic radius (jCantiello et al.ll2005l) . 

The targets analysed in this work have already been 
considered for other SBF measurements surveys. In par- 
ticular, the afore mentioned ground-based SBF survey 
(|Tonrv et al.l l200ll ) the near IR measurements for nine 
Virgo ellipticals by iPahre fc Mould! (1 19941) . and the ACS 
Virgo Cluster Survey (|C6te et al.l 12004) all have the two 
galaxies in their target list. The SBF magnitudes from these 
databases are reported in the first part of Table |3J None of 
these cited works presents the measurement of SBF magni- 
tudes at various galactic radii. On the contrary, thanks to 
the high quality of the VLT data available, we have been 
able to obtain fluctuation amplitudes at various radii (Table 

EJ. 

Differently fro m other bright elliptical gala xies studied 
in previous works (|Cantiello et alJl2005l l2007al) , we can re- 
port the absence of radial V and R gradients in NGC 4374. 
For NGC 4621, instead, a small but non-negligible bright- 
ening of SBF magnitudes at larger radii is found in V and 
R. On the other hand, the B of both galaxies shows a signif- 
icant scatter between the different annuli, thus no evidence 
of systematic trends with radius is observed in this band. 
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Fig. 4. SBF-colours and integrated-colour s meas ured for NGC 4621 (left panels) and NGC 4374 (middle panels), compared 
to the SPoT models from lRaimondo et al.l ([201 If ). Black symbols mark observational data for the various annuli as labelled 
(triangles for NGC 4621, squares for NGC 4374). Measures within regions hosting SNe, and the measure for the main 
annulus for both galaxies are shown in right panels: NGC 4621 and its SN la (SN 1939B) are shown with a full and empty 
triangle, respectively; NGC 4374 and its SNe la (SN1957B, and SN 1991bg) are shown with full and empty squares, 
respectively. For the models, different [Fe/H] are indicated with different symbols ([Fe/H] — -0.7, -0.4, 0.0,-1-0.4 dex 
with three, four, five and seven pointed stars, respectively). Larger symbols size mark older ages, from 4 to 14 Gyr with 
2 Gyr step. [See electronic version of the Journal for a colour version of the figure.] 



4.1. SBF-colours versus integrated colours 

By plotting SBF and integrated colours for each target 
(Figure |4]) , and comparing the data of the two galaxies be- 
tween them and with models, we can make various consid- 
erations on the host unresolved stellar syst ems. We consider 
the u pdated version of the SPoT models (jRaimondo et al.l 
1201 It ), which for the photometric bands and chemical com- 
position used in t his section confirm the results obtained 
from the previous iRaimondo et al.l (|2005f) models. 

As a first general consideration, we find that the SBF 
V— R colour, and the integrated B—V and V — R colours 
predicted by models representing old ( It 7 Gyr) and metal 
rich ([Fe/H] ^ —0.4 dex) SSPs are in good agreement with 
the measured values. Differently, the B is much more sen- 



sitive to the properties of the dominant stellar population 
(see below) and the matching with "canonical" SSP models 
is less satisfactory in this band. 

Starting from these general considerations on the prop- 
erties of the "dominant" stellar population in the galaxies 
(i.e. the stellar population which is is emitting the largest 
part for the flux responsible for the colours and the SBF 
signal) let us analyse more in details the result of Figure 0] 

The left and middle columns of panels in Figure HI show 
the SBF colours and V—R measurements in the five annuli 
considered for each galaxy versus the integrated B—V. The 
annuli appear numbered in the lower panels according to 
the numbers reported in Tabled so that inspecting the data 
in the figure it is also possible to recognise the radial be- 
haviour of plotted quantities. In the third column, instead, 



9 



Cantiello et al.: BVR SBF analysis for NGC4621 and NGC4374 



the overall SBF measures are plotted ( "main" annulus data 
in Table 02 empty symbols) together with SBF in the SNc 
regions. It is useful to emphasise that for both galaxies our 
integrated colours a nd colour profiles agree with the same 
measurements from lldiart etall (|2002t) . 

By inspecting results in the figure, we find that 
NGC462I does not show any significant SBF colour gra- 
dient either in B—V or V—R (lower and middle left pan- 
els), and no gradient even in the integrated V—R colour 
(top left panel). The only colour showing radial changes 
is the B—V, varying from ~ 0.91 mag for the inner an- 
nulus (#1), to ~ 0.84 in the outermost (#5). As already 
mentioned, inspecting the V and R data for NGC4621 in 
Table OS a small but non-negligible gradient seems to be 
present, of ~ 0.25 mag in V and _~ 0.15 mag in R, if the 
inner annulus is discarded. The V and R gradients have 
the same trend with colour, i.e. brighter SBF magnitudes 
associated with bluer colours/larger radii, thus they tend 
to cancel-out when the SBF-colour is considered. In fact, if 
the innermost annulus is neglected, a SBF-colour gradient 
of appears in the V—R versus B—V plane for NGC4621, 
although it is of the same amplitude of estimated uncer- 
tainties. Even so, the radial trends of the SBF magnitudes 
and of the B—V are both consistent with the well known 
scenario of a more metal-poor stellar component (i.e. bluer 
integrated colours, and brighter SBF magnitudes for the 
pass-bands considered here) at larger galactocentric radii. 

For NGC 4374 (central panels) we find that V-R and 
V—R are practically flat along the galaxy radius, if the 
innermost annulus (#1) is excluded. Due to the presence 
of a dusty disk in the galaxy centre, in fact, we left un- 
masked only a small fraction of the area in the annulus #1, 
but further contamination from dust might still be present. 
Differently from the previous case, the V and R data of 
NGC 4374 in Table [3]_do not show any systematicity versus 
radius, although the V—R versus B—V panel seems to show 
a correlation between the two plotted magnitudes, mostly 
due to the cited innermost annulus. 

As mentioned above, the comparison between measure- 
ments and models in the two upper panels of each column 
reveals that the data of galaxies lie close to the position of 
SSPs with [Fe/H] > - 0.4 dex and t > 7 Gyr. However the 
age-metallicity degeneracy suffered by classical and SBF 
colours obtained by combining these bands prevents us from 
giving more precise information on the dominant stellar 
components in the galaxy. 

The lower panels, which involve B, deserve a separate 
and detailed discussion, since both galaxies show a. B—V 
colour much bluer than models, an evidence more severe 
for NGC 4621. For this galaxy, in fact, the {B—V)Main is 
~0.2 mag bluer than NGC 4374, and ^0.4 mag bluer than 
models. We also note that NGC 4374 shows a radial change 
leading the B—V value of the outermost annulus to be 
quite similar to the models. Differently, NGC 4621 does not 
present any systematic gradient and, except for the inner- 
most annulus the B—V is nearly constant in every studied 
region of the galaxy. 

To understand the origin of the mismatch between data 
and models, visible in the lowermost panels of Figure[4j and 
of the peculiarly blue B—V colour of NGC 4621, we have 
considered to take advantage of the versatility of the SPoT 
stellar population synthesis code and carry out some spe- 
cific numerical simulations, to investigate if and how the 
presence of a complex stellar population modifies the ex- 



pected B—V . In addition, for metal rich models we have 
also considered the case of a non negligible fraction of 
Horizontal Branch (HB) stars having higher effective tem- 
peratures with respect to "canonical" stellar evolution mod- 
els. 

In order to find indications on the origin of the peculiar 
B—V behaviour of the two galaxies, our numerical tests 
have been organised as follows: 1) we have explored the ef- 
fects of the presence of a hot HB (HHB) component to field 
stars; 2) a young stellar population is added to a "reference" 
old and metal-rich population, and 3) a metal-poor popula- 
tion is added to the reference one. For the comparison in B, 
V and R bands we adopted the SBF measured in five annuli, 
plus the measures in the SNala regions, while for the panels 
including the K taken from literature we adopted our main 
annuli SBF measurements. The results of the simulations 
and the comparisons with data are briefly discusses below, 
and shown in Figure [5] where three different SBF colours 
are plotted against B—V . 

— In order to mimic the presence of HHB stars, we car- 
ried out a numerical experiment by considering that 
50% of HB stars have suffered large mass-loss dur- 
ing the RGB phase (first column of panels in Figure 
O. As a result, a percentage of low mass stars popu- 
lates the blue/hot part of the horizontal branch. This 
is not unrealistic since such hot component is observed 
in several metal-rich stellar clusters (e.g. NGC 6441 and 
NGC 6388) , whatever the mechanism responsible of this 
evidence is. In the figure we consider three populations 
with solar metallicity and ages 10, 12, and 14 Gyr. Both 
SBF and integrated colours move from the "standard" 
positions (standard HB, SHB, label in the figure) to the 
hot /blue region of the diagram (HHB in the figure) . The 
effect is larger for SBF colours including B with respect 
to integrated colour due to their increased sensitivity to 
even a small number of bright hot stars. 

— In the second numerical experiment (middle panels) a 
very young population (two different ages are consid- 
ered: 30Myr and 100 Myr) is added to a population of 
14 Gyr and solar metallicity, with mass ratio between 
the two components l id ■ 10yo Ung - From the middle 
panels of Figure [5] appears that a very recent and dif- 
fuse burst of star formation is required to obtain a good 
matching between models and data. 

— The third column of Figure [5] reports the effects of a 
secondary metal-poor component of 14 Gyr added to the 
main solar metallicity one, the [Fe/H] of the metal-poor 
component are labelled in the lower right panel of the 
Figure. Only a very metal poor population with a mass 
comparable to the one of the main component {l m .rich ■ 
lm.poor) produces relevant effects on SBF colours. In 
that case, however, the values of the observed integrated 
coluor B—V are not well reproduced any more. 

The three sets of simulations described seem to point 
out what already suggested: the most likely solution to the 
puzzling behaviour of £?-band SBF is the presence of hot 
HB stars. In fact, while the presence of a hot/young stellar 
component provides a good matching between integrated 
and SBF colours models with data (second scenario, middle 
panels in Figure [5]), it appears unrealistic that these regular 
ellipticals host such diffuse and very young stellar compo- 
nent. Moreover, in such a case, i.e. an object with massive 
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Fig. 5. SPoT stellar population models obtained with non-canonical assumptions. In all cases the initial (standard) 
population has solar metallicity and t~ 14 Gyr, and is shown with a blue star in each panel, while the final composite 
population is shown with a different colour and connected with a dashed line to the reference model. Upper panels show 
main annulus measures, middle and lower panels show multi-annuli measures and SNe la regions. Left panels: a fraction 
equal to 50% of the total HB-stars is simulated being HHB. For these models three ages are considered (t=10, 12 and 14 
Gyr, increasing ages are marked with larger symbols). Middle panels: a young population of t=30 Myr (red), or t=100 
Myr (magenta) is added to the old solar one. The fraction in mass of old to young stars is reported in the lower panel. 
Right panels: various metallicities are added to the solar one, as labelled. The mass fraction metal-poor to standard SSPs 
is shown in the lower panel. Symbols for observational data are as in Figure 2] [See electronic version of the Journal for 
a colour version of the figure.] 



and extended recent star formation, it is likely that a large 
quantity of dust would still be present, preventing the SBF 
measure itself. However, with the only exception of an in- 
ner r 20" dust ring in NGC 4374, we do not find any 
sign of extended dust in both galaxies. The two-metallicity 
mixing scenario seems even more unlikely, such mixing, in 
fact, does not solve the mismatch between data and models 
for optical bands. 

As a further element in support of the role played by 
HHB stars in deter mining the SB F signa l of NGC 4621, we 
recall the works by I Brown et al.l i|2l)l)l)i and iBrown et al.l 
(|2008l) , based on deep near and iax-UV images of the com- 
pact elliptical galaxy M32. Using HST data, the authors 



found that the number of PAGB stars in M32 is signifi- 
cantly lower than the expectations of their stellar evolution 
models, while the presence of a HHB population has been 
observed and identified as the main contributor to the UV- 

emission of the galaxy. 

Related to this issue, ICantiello etaH (j2007aft , using 
the B data of eight ellipticals, suggested that the role of 
hot evolved stars cannot be neglected in modelling SBF 
magnitudes in this pass— band. In that case,_though, the 
mismatch between the data and models in B was solved 
by enhancing the number of P AGB expected assu ming 
the evolutiona r y pres criptions bv lBrocato et al.l |l990) and 
iBrocato et all (feoOOft . 
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As far as concerns NGC4621, our present results point 
towards the direction suggested by the observations of M 32. 
We show that HHB stars can be the stellar component re- 
sponsible of the observed -B-excess, even though a combi- 
nation of two stellar components (HHB and PAGB stars) 
could not be ruled out and the topic requires further inves- 
tigations. 

As already mentioned for the case of M32, a further- 
piece of information comes from the comparison of the in- 
tegrated properties of the two galaxies in the wavelength 
interval where the very hot stellar component is dominant, 
i.e. th e UV regime. In this wavelength regime. lLongo et al.l 
(fl99l found that NGC4621 is brighter than NGC4374. 
The puzzling presence of a strong UV emission in some reg- 
ular early-type galaxies, d iscovered in late '70s (jCode et al.l 
U972tlBertola et all [1980 ). is now widely interpreted as the 
presence of a old hot stellar component. Although the mech- 
anisms that created such component , or its true nature, 
are not well und erstood (|Park fe Led Il997t iKavirai et al.l 
201)7: |Han| 120081) . some of these old hot stellar sources 
may have effects on B, as discussed by various authors 
based on both SSP models predictions (IWorthev 1993a ; 
Cantiello et al.l2003f). or on empirical data (IShopbell et all 



1993; ISodemann fe Thomsenlll996HCantiello et al.ll2007al). 

Moreo ver, a recent study of lBuzzoni fc Gonzalez-Lopezliral 
( 2008) presented a detailed analysis on the link between the 
UV characteristics and near-/i? SBF amplitudes of ellipti- 
cal galaxies. 

In conclusion, the present analysis seems to support a 
scenario where the peculiar B— V is related to a hot and old 
diffuse stellar component, like HHB stars. Larger samples 
of SBF colours in blue bands are required to provide further 
constraints. 



4.2. SBF colour-colour diagrams 

The coupling of literature /, z and X-band SBF data for 
the two galaxies with our measures (main annuli) provides 
a sample of six independent SBF values useful to obtain 
three independent SBF-colour versus SBF-colour diagrams. 
Such wide pass-band coverage is, to our knowledge, the first 
ever presented for SBF. The diagrams with observations 
are shown in Figure [51 together with models. The upper 
panel of the F igure shows a fe ature that has alread y bee n 
discussed (e.g. iBlakeslee et al.|[200ll: ICantiello et al.ll2003h , 
that is that model separation for pure optical SBF colours 
is not effective in removing the age-metallicity degeneracy. 
In such panel, in fact, all [Fe/H] >-0.4 dex models overlap 
to each other, with old and metal rich models superposed 
to young and metal poor ones. 

The models degeneracy is significantly reduced in the 
middle panel, which includes the optical to near-IR V—K 
SBF colour. The position of the two objects overlaps with 
models with different metal contents. In particular the 
bluer point, i.e. NGC4621, is located near to the region 
of [Fe/i2]~-0.4 dex with respect to NGC4374, which ap- 
pears slightly more metal-rich. Such behaviour is not un- 
expected, since the latter galaxy is brighter than the first 
one and, due to the known mass-metallicity relation (e.g. 
iBernardi et al.l[l998f) , it is reasonable to expect a field pop- 
ulation slightly more metal-rich in NGC 4374. It is worth 
to emphasise that the two- SBF-colour panel represents a 
useful tool to provide insights to the "absolute" [Fe/H] of 
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Fig. 6. SBF colour-colour diagrams obtained from present 
measures (main annulus) and literature data. Symbols and 
models are the same as for Figure [4] For consistency 
with other filters the z-band SBF magnitudes given by 
Bla keslee et al.l ((2009) are transformed from the ABmag 
to the Vegamag system according to lSirianni et al.l ((2005). 
[See electronic version of the Journal for a colour version 
of the figure.] 
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the dominant stellar component, within the observational 
error-bars, and the SSP models scenario adopted. 

In other words, based on the SBF-colours analysis and 
on t he integrated colours of field stars from this work (see 
also iTonrv etaUfcOOlt lldiart et alJ l2002t iMei et alJl2007h 
we find that the dominant stellar population of NGC 4374 
is either more metal rich or older, or both, than what found 
in NGC 4621. Moreover, this result agrees with independent 
estimates of the metal content in these galaxies, based on 
spect ral index measure ments, which predict nearly iden- 
tical (lldiart et al.1 1200^) or slightly higher metallici ty for 
NGC 4621 (|Kuntschner et al.ll200lfc IScott et al.ll2009h . 

It must be pointed out, though, that these results are 
based on the naive interpretation of field stellar popula- 
tion properties solely in terms of age and chemical content 
differences. Nevertheless, as discussed above, SBF magni- 
tudes and colours can be very sensitive to the presence of a 
blue/hot stellar sub-system, especially SBF in bluer pass- 
bands. 

The lower panel in Figure [6] includes our B measure- 
ments. As for the first panel, the V—z colour is used in 
abscissa, in this case the separation between different SSP 
models allows to recognise that for the two larger [Fe/H] 
values the V—z colour is bluer for older ages. Taking as 
reference the models with solar metallicity, one can con- 
clude that both galaxies have similar V—z, matching with 
the oldest SSP models, but NGC 4621 has B—K ~ 0.6 mag 
bluer than NGC 4374, as one would also expect in the case 
of a diffuse hot stellar component is present in this galaxy. 
As an example, a canonical SSP model at t — 14 Gyr, solar 
metallicity has B—K ~ 8.4 and V—z ~ 3.7 mag; by artifi- 
cially increasing the number of HHB stars (as described in 
section B~Tj) our models predict B—K ~ 7.7, i.e. roughly the 
AB—K colour difference between the two galaxies, while 
V—z is left unchanged. 

A further consideration is the fact that in the bottom 
panel of Figure[Bl NGC 4621 falls within the region of mod- 
els with [Fe/H] between -0.7 and -0.4 dex, even though 
the £?-band SBF is too bright. Such behaviour could be 
explained with bright i-T-band SBF, so that the B and K 
SBF excess compensate each other allowing observations to 
match with models of intermediate metallicity. This would 
also make V — K too red, causing the shift observed in 
the middle panel. In that case, i.e. too bright B and K- 
band SBF, the three panels of Figure [6] are all suggesting a 
[Fe/H] < - 0.4 for NGC 4621, and a larger metallicity for 
the dominant stellar component in NGC 4374. 

Whether the position of NGC 4621 with respect to mod- 
els is due to the presence of a metal poor dominant stel- 
lar component, or it is related to a diffuse hot component, 
cannot be established with present data. New independent 
observational data sets, and accurate modelling, able to 
match all the observed properties of the two galaxies, are 
needed to better understand the physical characteristics of 
the unresolved stellar systems in both targets. 

Nevertheless, the coupling of £?-band with other opti- 
cal or near-/i? SBF measures, appears to be a promising 
method to unveil the properties of hot stellar components 
possibly hidden to other photometric indicators. 



4.3. SBF colours and stellar populations properties in regions 
hosting type la Supernovae 

The correlation between the SN la peak magnitude and the 
host-galaxy morphology, i.e. the host stellar pop ulation, is 
well known since the work bv lHamuv et al.l (Qjj96;), showing 
that intrinsically fainter events occur in early-type galaxies, 
while luminous events are often seen in late-type galaxies. 
The results on the mean stellar population properties in the 
galaxies, discussed in the previous sections, seem to support 
the idea that old stellar system, i.e. old progenitors, are 
required for sub-luminous SNe la. 

To better constrain the properties of the local stellar 
populations we have measured SBF and colours in the areas 
where the SNe la exploded, adopting a 25" x 25" square 
region. The results of the measurements are reported in 
Table[3j and shown with full triangles (for NGC 4621), and 
full squares (for NGC 4374) in the right panels of Figure |3J 

First of all, we note that, from B—V, V—R and V—R 
data, the stellar population properties in the selected re- 
gions and in the main annulus (empty black symbols) ap- 
pear remarkably homogeneous. Taking into account only 
such colours, we conclude that the stellar component in the 
regions hosting the three SNe la and the overall stellar pop- 
ulation in the galaxies have quite similar properties, with 
all data matching SSP models of [Fe/H] <i — 0.4 dex and 
t ^ 7 Gyr. In any case, old ages are ex pected, support- 
ing the results by iGallaeher et al.1 (|2008h based on spec- 
troscopic data. It is interesting to note that the properties 
of the dominant stellar populations in the regions of the 
two SNe la host by NGC 4374 appear strikingly similar in 
these panels, nearly identical to each other, notwithstand- 
ing the large projected separation (~10 kpc) between the 
two regions. 

Larger differences are seen between the three regions in 
the B—V versus B—V panel (lower right panel in Figure 
E]). Due to the aforementioned stronger sensitivity of B- 
band SBF to changes of stellar population properties re- 
spect to other pass-bands, the larger scatter between the 
three SNe la is not surprising, and it might be possibly 
related to different levels of local "pollution" from a hot 
stellar component. Furthermore, the regions of NG_C4374 
where SN 1991bg and SN 1957B exploded show a AB—V ~ 
0.4 mag,_despite their nearly homogeneous B—V, V—R 
and V—R colours. Whether the AB—V is related to the 
SNe progenitors and their environment, and thus to the 
AB. max ~ 2.3 mag between the two SNe la, or to problems 
with the P r correction cannot be said using present data, 
and further analysis based on a richer sample is needed. 

For what concerns NGC 4621, even including B-band 
SBF data we find a good matching between global and lo- 
cal, i.e. near SN 1939B, stellar population properties. 

5. Conclusions 

We have presented a detailed multi-band study of SBF 
magnitudes for two bright galaxies in the Virgo cluster, 
NGC 4621 and NGC 4374, based on deep B, V and i?-band 
imaging data of the FORS2 camera at the VLT telescope. 

Among the three bands available, the V and R have 
known empirical SBF calibrations useful to obtain the dis- 
tances of the targets. Coupling our measurements with such 
calibrations, or with calibrations obtained from simple and 
composite stellar population models, we obtained accurate 
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distances for the two galaxies which agree very well with 
other distance estimates taken from literature. This demon- 
strates both the reliability of the measured SBF magnitudes 
and the goodness of the calibrations adopted. Taking ad- 
vantage of the fact that various SBF distances, based on dif- 
ferent calibrations for different filters, are available for the 
two targets, we have compared the average SBF distance 
moduli with other non-SBF ones, to check for any possible 
systematics. Although the sample of galaxies is statistically 
small, the results obtained seem to rule out the presence of 
possible bias on SBF distances. 

Taking the median of all available SBF and non-SBF 
distance estimates, we estimated: m — M = 30.98 ±0.17 for 
NGC4621, and m-M = 31.22 ± 0.16 for NGC4374. 

The two target galaxies hosted a total of three SNIa 
events: one in NGC 4621, SN 1939B, and two in NGC 4374, 
SN 1957B and SN1991bg, all classified as sub-luminous. 
Using some recent calibrati on relations based o n decline 
rate and colours of the SN (iFolatelli et al. 12010ft . and on 
the updated MLCS method (|Jha et al.ll2007l ) we find an ex- 
cellent agreement between the SBF and SNela distances. 
A result that is very promising in view of a reduction of 
the number of "rungs" to bridge local to cosmological scale 
distances, i.e. to significantly reduce the systematic uncer- 
tainties on distances of objects at large distances and, con- 
sequently, on cosmological parameters. 

We also carried out SBF measurements on _B-band 
images, but such measures were rejected for distance 
analysis. The sensitivity of SBF magnitudes in this band 
to the stellar population properties makes unreliable any 
calibration in this band. We tentatively derived a B-band 
calibration based on SPoT models - which have proved 
being realistic i n both V and R ban ds, as well as other 
studied bands (jBiscardi et al.l 12008ft -, but the results 
obtained confirm that SBF magnitudes in such band are 
not well suited for any kind of distance analysis. 

Thanks to the high quality of the data, we have suc- 
ceeded in measuring SBF magnitudes in various galactic 
regions. Both integrated and SBF colours (excluding B) 
seem to point out that the stellar population are relatively 
uniform along the galaxy radius. No sizable SBF gradient 
is observed in NGC 4374, while for NGC 4621 a small but 
non-negligible SBF gradient in V and R is observed, con- 
sistent with the known scheme of more metal poor domi- 
nant stellar populati ons at larger galactocentric radii (e.g. 
ICantiello et al.ll2005l ). Given the fact that i) compared to I- 
band SBF V and R are more sensitive to stellar popul ation 
properties (|Blakeslee et al.ll200 It ICantiello et al.1 12003) . and 
ii) I gradients up to ~ 0.4 magnitudes were measured by 
ICantiello et al.l (|2005ft on similar radial scales for a differ- 
ent set of targets, the negligible V and R gradients detected 
here might be explained by the position of the two objects 
within their host cluster. Since the two targets analysed 
are located in the regions with highest galaxy density in the 
Virgo cluster, in fact, this behaviour might be related to the 
major-merging events in such environment, wh ich tend to 
smea r out possible stellar population gradients ()Kobavashil 

Taking advantage of the multi-wavelength coverage of 
our data set, and of existing /, z, and if-band SBF mea- 
surements, we analysed the data of both galaxies in the 
SBF-colour versus integrated-colour and SBF-colour ver- 
sus SBF-colour diagrams. This is the first SBF-analysis 



carried out using SBF data in six different pass-bands. As 
a result we concluded that the dominant stellar component 
in the two galaxies is very similar, though NGC 4374 seems 
to be slightly more metal rich than NGC 4621. 

If B-band measures are taken into account, the SBF 
colours of the two galaxies show non- negligible differences, 
with NGC 4621 having brighter SBF than NGC 4374. Given 
the known link between SBF magnitudes in blue b ands and 
the properties of a hot old stellar component ()Worthevl 
Il993at ICantiello et "all 12003ft . we used the SPoT Stellar 
Population Synthesis code to simulate populations with 
"non-canonical" properties. In particular, to a solar metal- 
licity old t ~ 14 Gyr population, we have 1) enhanced the 
content of hot HB stars (i.e., ~ 50% canonical HB and ~ 
50% HHB), 2) added a very young (down to 30 Myr) dif- 
fuse secondary component, and 3) added a more metal poor 

SSP (down to [Fe/H] 2.3 dex). Within the limits of the 

small number of data, and adopting as a razor the physical 
plausibility of the population mixing, the simulations seem 
to favour the HHB component scenario. 

Th is is also support e d by the fact that iBrown et al.l 
(|2000D and IBrown et all (|2008ft find the presence of HHB 
stars in the compact elliptical galaxy M32 on the basis 
of UV observations. In this framework, it is relevant to re- 
call that NGC 4621 is substantia l ly bri ghter than NGC 4374 
in the UV bands (|Longo et al.l 1 199 ll ). and that the inte- 
grated properties of normal elliptical galaxies with bright 
UV emission are now int erpreted as an effe ct of a diffuse 
hot stellar component fe.g. lPark k, Ledfl997ft . Thus, further 
analysis is required to better constrain the links between the 
UV emission and SBF in early type galaxies, also in view 
of the recent stu dies on the relation of UV with ne& i-IR 
SBF amplitudes ()Buzzoni fc Gonzalez-L6pezliral20 08) , and 
larger samples of SBF colours in blue bands are required to 
provide further constraints to this scenario. 

Finally, taking into account the integrated colours and 
V plus R SBF for the areas of the galaxies where type la 
Supernova exploded, we find no substantial differences be- 
tween the local and the global stellar population properties 
in the galaxies. For SN1939B and its host, NGC 4621, this 
is also true for JS-band SBF data. On the contrary, there 
is a AB—V ~0.4 mag difference between the regions of 
the two SNe host in NGC 4374, SN 1957B and SN1991bg. 
Due to the quoted sensitivity of blue-band SBF to stellar 
population properties, such difference and, more in general, 
the scatter in B— V between the three SNela might be re- 
lated to different levels of local pollution from a hot stellar 
component, although for the case of SN 1991bg we cannot 
exclude the possible systematic effect of a low B signal, 
comparable to the P r correction. 

In conclusion, the present set of SBF measures shows 
that optical to near-/i? SBF magnitudes can be very effec- 
tive to unveil the properties of global and/or local stellar 
populations in distant galaxies. In particular, if B (or even 
SBF in bluer bands) are available, various constraints can 
be set to the role of hot field stars in normal ellipticals, 
at the same time multi-band SBF can be used to anal- 
yse the relation between SNe and local stellar population 
properties. Nevertheless, serious limitations are set by the 
small samples of measures still available, new observations 
are advisable to further study stellar population properties 
based on multi-colour SBF analysis. 
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